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Abstract 

With the establishment of the AUT University 12m radio telescope at Warkworth, New Zealand has now 
become a part of the international Very Long Baseline Interferometry (VLBI) community. A major 
product of VLBI observations are images in the radio domain of astronomical objects such as Active 
Galactic Nuclei (AGN). Using large geographical separations between radio antennas, very high angular 
resolution can be achieved. Detailed images can be created using the technique of VLBI Earth Rotation 
Aperture Synthesis. We review the current process of VLBI radio imaging. In addition we model VLBI 
configurations using the Warkworth telescope, AuScope (a new array of three 12m antennas in Australia) 
and the Australian Square Kilometre Array Pathfinder (ASKAP) array currently under construction in 
Western Australia, and discuss how the configuration of these arrays affects the quality of images. Recent 
imaging results that demonstrate the modeled improvements from inclusion of the AUT and first ASKAP 
telescope in the Australian Long Baseline Array (LB A) are presented. 
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1 Introduction 

The Warkworth Radio Astronomical Observatory 
(WRAO) is located some 60 km north of the city of 
Auckland, near the township of Warkworth. The 
observatory is operated by the Institute for Radio 
Astronomy and Space Research (IRASR) of AUT 
University. The observatory's 12-m radio telescope 
operates in three frequency bands centred on 1.4, 
2.3 and 8.6 GHz. In addition to astrophysical 
observations this fast-slewing ( 5 ° per second in 
Azimuth) antenna is well suited to the purposes of 
geodetic VLBI [1^ and spacecraft navigation and 
tracking [2]. 

In February 2011 the AUT University 12m ra- 
dio antenna officially joined the Australian Long 
Baseline Array (LBA) and now regularly particip- 
ates in its VLBI sessions. The primary product 
of this VLBI work is high resolution radio domain 
images from which the physical properties of radio 
astronomical sources are studied. 

The AUT University 12m radio antenna has ex- 
panded the maximum east-west baselines of the 
LBA [3^ by almost a factor of two, from ^ 1300km 
to ^ 2400km. With the further addition of the 
ASKAP antennas in Western Australia east-west 
baselines of ^ 5000 km are achieved, providing 
corresponding increases in resolution. The goal 



of this paper is to investigate the way for further 
improvement of the quality of the array, first of all 
in terms of extension of the north-south baselines 
to yet further improve the resultant image quality. 

In Section 2, we briefly present the theory behind 
radio interferometry and review the current meth- 
ods used to obtain radio images of astronomical 
sources. Section 3 specifically outlines current 
image recovery methods. In Section 4, we model 
the effects of adding the AUT University 12m radio 
antenna which will be referred too as Warkworth, 
the new ASKAP [4 and AuScope antennas 
on the imaging performance of the LBA. Section 
5 discusses current activities and presents actual 
images achieved with the AUT 12m antenna as 
part of the LBA. 

2 Radio interferometry 

A radio interferometer is a pair (or more) of 
antennas used to measure the visibility function 
due to the sky brightness within the field of view 
of the antennas [6j. From sampling this visibility 
function it is possible to recover an image of the 
observed field of view [6^. 



2.1 UV coverage for models 

Imaging of a radio astronomical source by the tech- 
nique of aperture synthesis was first demonstrated 
by Prof. Ryle [7 using the Cambridge Radio 
Telescope [8^. Antennas able to track a source for 
an extended period as the Earth rotates will trace 
out elliptical paths in the u^v plane (orthogonal 
plane to the direction of the astronomical source). 
The components u^v^w may be determined from 
the expression [6 : 
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Here 5o are the hour angle and declination of the 
source, A is the wavelength of the radio frequency 
being observed. It is customary in VLBI to 
eliminate the Hq (hour angle) term by setting the 
X axis of the coordinate system to the Greenwich 
meridian Hq = 0, resulting in: 
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As Lx^Ly and are constants for a given pair of 
antennas, this is the equation of an ellipse in the 
u,v plane (it becomes the equation of a circle u^-\-v'^ 
when So = 90°). For an array of N antennas we 
will have N{N — l)/2 pairs of elliptical loci. 

Plots of these loci (tracks) demonstrate the pro- 
gressive improvement in filling of the v plane as 
additional antennas are added to the LBA array: 
Warkworth plus LBA in Figure [l] Warkworth, 
LBA and ASKAP in Figure [2) Warkworth, LBA, 
ASKAP and AuScope in Figure |3] All plots are 
for a common source declination of —37°, that is 
for a Southern Celestial Hemisphere radio source. 
Both u and v are given in wavelengths, A. 
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Figure 1: UV Coverage for Warkworth and the LBA 

2.2 Visibility function 

The astronomical source is treated as a two dimen- 
sional image of intensity /(/,m) on the celestial 
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Figure 2: UV Coverage for Warkworth, LBA and 
ASKAP 



ww-lba-askap-aus at 1.420 GHz in XX 2011 Aug 09 
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Figure 3: UV Coverage for Warkworth, LBA, ASKAP 
and AuScope 

sphere, / and m being coordinates in the plane of 
the celestial sphere. The projection of / and m 
to the plane perpendicular to the direction of the 
astronomical source from the Earth, as defined by 
the coordinate system u and is shown in Figure 

m 

It can thus be seen that the visibility measured 
by an interferometer is a sample of the visibility 
function V{u^v), and may be expressed as the 
Fourier transform of the modified sky intensity / 
i: 



V{u,v) = j j I{l,m)e-^''^^''^^'''^Uldm (2) 

As the visibility is sampled out to a maximum 
radius bmax in the v plane, the array produces 
information similar to a single circular aperture of 
diameter D ^ ^bmax- Resolution of the image 
will be ^ l/bmax^ but the quality of the image 
is determined by the v coverage. 



m 




Figure 4: Geometric relationship between the source 
observed /(/, m) and an interferometer or antenna pair 
upon the u, v plane [9 

3 Image recovery methods 

During an observation, sources are tracked by 
the interferometer array and the signal at each 
telescope separately recorded. Data from all 
telescopes is then sent to a processing centre for 
correlation. As F is a Fourier transform of the 
source brightness distribution / Equation ([2|, the 
latter can be recovered by means of the inverse 
Fourier transformation: 

/(/,m) = J Jv{u,v)e^'''^''^^'''^Uudv (3) 

As integration occurs over u and the more 
information we have about V{u^v)^ the more fully 
we can conduct numerical integration in[3j In other 
words: the better the u^v coverage, the greater 
chance to restore the image, I(lm)^ without errors 
and artefacts. 

Prior to the inverse Fourier transforming [3j calib- 
ration steps are invariably required as variations in 
performance and characteristics of antennas must 
be accounted for (the array is not homogeneous). 
Additionally, data may have to be flagged due to 
factors such as equipment failure at an antenna or 
occurrence of radio frequency interference (RFI). 

To take advantage of the computational efficiency 
of the Fast Fourier Transform (FFT), the correl- 
ated visibilities must be binned onto a regular 
grid. Weighting schemes may be applied at this 



point that influence this gridding and aspects of 
the fidelity of the final image [lOj. Once "gridded" 
the data is Fourier transformed to obtain a first cut 
image. This first image is commonly referred to in 
radio astronomy as the "dirty map" as it consists 
of the true image of the source convolved with 
(corrupted by) the point spread function (PSF) 
of the imaging instrument, the VLBI array itself. 
Through use of one of various algorithms noted in 
the following subsection, an approximation to the 
true image is sought through deconvolution of this 
dirty map. 

3.1 Algorithms 

It was shown in Section 2 that the u,v plane is 
incompletely sampled by the antennas, thus an 
analytic approach cannot be used in recovering 
the true image. Several numerical algorithms have 
been found to effect a solution and are commonly 
used in radio astronomy interferometry to recover 
an image; CLEAN, MEM, NNLS. We only briefly 
present these different algorithms (more complete 
explanation and comparison is available from 
Cornwell [TO^) with more detail about CLEAN the 
most commonly used in VLBI imaging. 

CLEAN The most widely used radio astronom- 
ical image deconvolution technique is the CLEAN 
algorithm. Devised by H0gbom [11], this is a 
numerical deconvolution process applied in the 
(/,m) domain. The basic assumption is that the 
true image may be modeled as a set of point 
sources. In an iterative process the amplitudes 
and positions of these point sources are established 
through fitting the PSF of the VLBI array to the 
brightest point on the image. A final image is 
produced by convolution of this point source model 
with an appropriately sized Gaussian restoring 
PSF. There have been more recent refinements to 
this algorithm such as multi-scale CLEAN \12\ . 

MEM The maximum entropy deconvolution al- 
gorithm [T3^. 

NNLS Non-Negative Least Squares (NNLS) is 
discussed in some detail by Briggs [14 . 

3.2 Future trends 

Compressive sampling Until recently it has 
been held that to recover a signal or image the 
sampling rate must be at least twice the maximum 
frequency present in the signal [15 , at or above the 
Nyquist rate. A result of the newly developing field 
of compressive sampling finds that signals sampled 
at less than the Nyquist rate can be recovered if the 
information is "sparse". It was shown in Section 
2 that the v plane is not completely sampled 
by the antennas and thus we have a sparse model. 
Recently, some work has been undertaken to see 



if Compressive sampling can be used to recover 
images in the radio domain [16 . Feng Li et al 
[16] are attempting to implement these methods for 
deconvolving images in semi-real time. If this could 
be conducted in sequence with eVLBI correlation, 
then for the first time images could be available in 
near real time after an observation. 

3.3 Packages 

There are several software packages widely used 
by the Radio Astronomy community incorporating 
one or more of the above algorithms. These 
are the Astronomical Image Processing System 
(AIPS) Hi], MIRIAD [18^ and the standalone 
package DIFMAP for producing images from radio 
interferometers using the CLEAN deconvolution 
[19] algorithm. 

For the next generation of radio astronomical 
telescopes the Common Astronomy Software Ap- 
plications (CAS A) [20] is being developed by an in- 
ternational consortium, built on C++ tools under 
an iPython interface. This includes in addition 
to CLEAN and MEM a module for deconvolution 
using NNLS. 

4 Modeling NZ involvement in 
VLB! 

4.1 Modeling 

Consider the model image in Figure [5] of three 
elliptical Gaussian sources to represent a radio 
galaxy triplet. 
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Figure 5: Galaxy triplet model image 

This model image was input to a process developed 
at AUT[2T| that simulates and models VLBI arrays 
using a combination of the packages MIRIAD and 
DIFMAP. First taking the array configuration of 
Warkworth and the LBA we obtain the image 
presented in Figure (6] Without prior knowledge 
of the model image we would most likely conclude 



that we had a double source with one large ex- 
tended component to the right, a second more 
compact component to the left. Note that the 
restoring beam shape (bottom left hand corner 
solid shaded ellipse) infiuences the final image 
obtained. 

Clean XX map. Array; H0BART26M 
ww-lba.ant.uv at 1.420 GHz 2011 Aug 16 
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Figure 6: CLEAN Image from Warkworth LBA 

Following inclusion of an ASKAP antenna into 
the simulation as shown in Figure [7| it is just 
possible to resolve three separate components of 
the source, although the separation between two 
of the components is minimal. The beam (bottom 
left corner of image) is clearly more compact as a 
result of the additional resolution available. 

Clean XX map. Array: H0BART26M 
ww-lba-asl<ap.ant at 1.420 GHz 2011 Aug 16 
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Figure 7: CLEAN Image from Warkworth, LBA and 
ASKAP 

Finally, inclusion of the AuScope antennas adds 
additional and longer n-s baselines. A significantly 
more compact, although not quite circular beam 
results and three unambiguously separate sources 
are resolved as shown in Figure [8) 



Figure 8: CLEAN Image from Warkworth, LBA, 
ASKAP and AuScope 

5 Recent results 
5.1 VLBI 

The AUT University 12m radio antenna had been 
involved with the LBA prior to 2011 on an ad- hoc 
basis, but is now routinely included in experiments. 
First LBA science results that include the AUT 
and ASKAP antennas were achieved on April 2010 
with images of PKS 1934-638 obtained at 1.4 GHz 
[22] as shown in Figure [9j For comparison to 
demonstrate the improvement in imaging achieved 
compare this with an image taken of the same 
source using the LBA configuration before the 



addition of Warkworth and ASKAP Figure 10 




Figure 10: VLBI Image of 1934-638 at 1.4 GHz by the 

LBA, Credit: S. Tingay 

5.2 eVLBI 

eVLBI involves sending the received data from all 
antennas in real-time over a high speed network 
to the correlator. In July 2011 the first such LBA 
observations that included the AUT 12m antenna 
and newly commissioned ASKAP antenna took 
place with a source observed in L-Band (1.4 GHz), 
with the resultant image shown in Figure [TT] The 
data from the AUT University 12m radio antenna 
was streamed across the Kiwi Advanced Research 
and Education Network (KAREN) in real time at 
a rate of 512 Mbps to the correlation centre at 
Curt in University, WA. 




Figure 11: eVLBI Image of 0637-752 at 1.4 GHz by 
Warkworth, LBA and ASKAP July 2011, Credit: S. 
Tingay 



Figure 9: VLBI Image of 1934-638 at 1.4 GHz using 
LBA, Warkworth and ASKAP, Credit: S. Tingay 
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